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We revisit recent claims about the instability of non-rotating tunnel coupled annular Bose-Einstein conden- 
sates leading to the emergence of angular-momentum Josephson oscillation [Phys. Rev. Lett. 98, 050401 
(2007)]. It was predicted that all stationary states with uniform density become unstable in certain parameter 
regimes. By careful analysis, we arrive at a different conclusion. We show that there is a stable non-rotating 
and uniform ground state for any value of the tunnel coupling and repulsive interactions. The instability of 
an excited state with it phase difference between the condensates can be interpreted in terms of the familiar 
snake instability. We further discuss the sign of the tunnel coupling through a separating barrier, which carries 
significance for the nature of the stationary states. It is found to always be negative for physical reasons. 

PACS numbers: 03.75.Lm, 67.85.-d 



I. INTRODUCTION 

Bose-Einstein condensates (BECs) located in different min- 
ima of an external potential created by magnetic or light forces 
and coupled by tunneling through a potential barrier have 
been the host to many exciting developments and discover- 
ies in recent years |Q]]. Phenomena explored include analogs 
of the Josephson effect in double or multiple quantum-well 
structures ||2|, |3[ |4(], gap solitons of repulsive BECs Q5|] , and 
quantum phase transitions [6]. Often these systems are mod- 
eled by considering just one mode per potential minimum and 
their linear coupling provided by tunnling through a separat- 
ing barrier. These simplified models, which are usually la- 
belled as two-mode or multiple-mode models, variants of the 
Bose-Hubbard model, or the discrete nonlinear Schrodinger 
equation, are tailored to describe certain properties or aspects 
of the dynamics of the many-body system under investiga- 
tion. Although there is an abundance of literature on such 
models HSilEIIlllIlHQlEIIl, there still ap- 
pear to exist misconceptions about the nature of the effective 
model parameters, especially the sign of the tunnel coupling, 
as only few works attempt to calculate such parameters based 
on a more complete theoretical treatment fill fl3i [l4ll . 

The sign of the tunnel coupling bears special significance 
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in systems where the tunneling appears over an extended 
(at least one-dimensional) region of space. Such systems 
have recently been analysed by Bouchoule |[l7jl and Kaurov 
and Kuklov 01811 . who studied two parallel tunnel-coupled 
cigar-shaped BECs. In another recent work, Lesanovsky and 
von Klitzing studied the stability of tunnel-coupled annular 
BECs 1 19]. The latter paper points to an interesting dynamical 
instability leading to the spontaneous formation of angular- 
momentum fluctuations. We will show further below that the 
sign on the tunnel coupling bears consequences on the nature 
and stability of the stationary states found in the mean-field 
treatment of tunnel-coupled BECs. Specifically, we find that 
the system studied by Lesanovsky and von Klitzing has a sta- 
ble ground state for any value of the tunnel coupling and re- 
pulsive interactions. The instability of an excited state with it 
phase difference between the condensates can be interpreted 
in terms of the familiar snake instability l2(ill . The ground 
state of a rotating co-planar double-ring system is discussed 
inRef. HIl]. 



We examine the stationary states of double-ring BECs in 
Sec. [n] A careful analysis of the sign of the tunnel coupling 
used in effective models for BECs in double-well traps fol- 
lows in Sec. [Ill] Conclusions are presented in Sec. [IV] 
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II. STABILITY OF STATIONARY STATES IN 
DOUBLE-RING BECS 

The calculation performed in Ref. lll9ll starts from a num- 
ber of generally reasonable assumptions. Under the condition 
that radial excitations of the vertically stacked annular BECs 
are suppressed by the trapping potentials and the only mech- 
anism for coupling the two systems is via tunneling through 
a potential barrier, the Gross-Pitaevskii equation for the two- 
mode spinor wave function (x„, Xd) specialises to 



idrXu/d = -dlxu/d ~ \k\ Xd/u + 7 IXu/d| X' 



u/d 



(1) 



Here Xu(&) is the condensate wave function for atoms in the 
upper (lower) ring. Our equation (fl3 agrees with Eq. (2) of 
Ref. 01911 . except that we explicitly indicate the negative sign 
of the tunnel coupling. We give detailed reasons for the rel- 
evance of the sign of the tunnel coupling below in Sec. [TTTJ 
where we also show that the tunnel coupling is indeed neg- 
ative. At this point, we only note that the tunnel coupling k 
was assumed to be positive in Ref. [19] (see their Fig. 1), in 
contradiction to our findings. 

The most general form of the polar-angle-dependent 
wave function can be written as a Fourier series, Xu/d = 
(2n)~i Y^ m a m d) e lmip . Inserting this Ansatz into Eq. © and 
equating coefficients of the orthogonal Fourier components, 
we find 

id T of m ' = m OL~ m ' — |k| a^i ^+77^; ^""^ o4 ' a n ' Q ; m _ n+n /. 

This result differs from the corresponding Eq. (3) in Ref. lfl9^ 
in the tunnel coupling and the non-linear term. 

As a first approximation, it is reasonable to assume that 
only the m = mode is occupied in each of the two annuli. 
Straightforward calculation yields the new ground and excited 
state of the coupled-annuli system, which are the symmetric 
and antisymmetric superpositions of single-well states having 
chemical potential fj,± = £=p |«|, respectively, s = jNo/(2ir) 
is defined in terms of the equal number of atoms No in each 
well as in Ref. [19]. In order to study the stability of these 
states, finite but small amplitudes in the m ^ modes are 
assumed: 



(U/d) _ -V±T 



.(u/d) , 



(3) 



Here the subscript ± distinguishes perturbations to the ground 
and excited states, respectively. Inserting the perturbation (0 
into Eq. (O and linearising in the small amplitudes u, v yields 



(u/d) 
«m,± = 

(u/d) 

y j- = 



■E±\K\) 

■e±\k\) 



,(u/d) , 

(u/d) 
>-m,± 



■ £v -m,± ~ \ K \ u m,± > 
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(u/d) 



(d/u) 
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(4) 



The upper (lower) sign refers to the symmetric ground (an- 
tisymmetric excited) state. Crucial differences between our 
Eq. © and Eq. (5) in Ref. lfl9ll result in markedly different 



excitation spectra. We find that both the symmetric (ground) 
state and antisymmetric (excited) state share one branch, 



w\ = \J (m 2 + eY 



(5a) 



whose frequency is independent of the tunnel coupling. This 
was also found in Ref. 11911 . In contrast to these authors, how- 
ever, we find that the second branch differs for the two states: 



W2,± = V ( m2 + £ ± 2|k|) — e 2 



(5b) 



Clearly, ll>2,+ is always real for repulsive BECs (e > 0), im- 
plying stability of the symmetric (ground) state of the coupled 
annular condensates. In contrast, the antisymmetric (excited) 
state will become unstable for e > \n\ — m 2 /2 > 0, sig- 
nified by u>2,- becoming imaginary in this range. Our own 
numerical simulations of the time evolution of the antisym- 
metric state seeded with a small amount of noise show the de- 
velopment of angular- momentum Josephson junctions similar 
to those shown in Figs. 2 and 3 of Ref. lfl9ll . 

In attractive condensates where e < 0, imaginary solutions 
of u)x for 2e < — rn 2 indicate the well-known modulational in- 
stability towards the formation of localized peaks (bright soli- 
tons) in the individual rings. For the symmetric state, LU2.+ 
does not add new instabilities (with imaginary solutions for 
e < —m 2 /2 — \k\). The antisymmetric state, however, is 
further destabilised by the tunnel coupling due to imaginary 
frequencies of u>2- at e < \n\ — m 2 /2 < 0. 



III. SIGN OF THE TUNNEL COUPLING 

In our analysis so far we have assumed that the sign of the 
coupling constant n is negative. This lead to the symmetric 



state with aQ d) 



„(u) 



const • e v + T and a^ 



with 



fi + = e — \k\ being the ground state. Let us now briefly con- 
sider the consequences of the (hypothetical) case of a positive 
coupling constant k > 0. The analysis of SecUJcanbe carried 
out the same way as before, with the difference that \ k\ should 
be replaced by — |k| in all formulae. It is easily seen that, in 

will be the 



(d) 



-a, 



(u) 



this case, the antisymmetric state with a { 
ground state. Since the sign change also affects Eq. ( Bbl ). we 
find the antisymmetric state being stable (for e > 0) and the 
symmetric one becoming unstable. However, since the roles 
of these states have changed, we still find that the ground state 
is stable for repulsive BECs. 

In order to determine the correct sign and value of the 
coupling constant k appearing in Eq. (Q]), we briefly revisit 
the derivation of this model. For the purpose of finding k, 
the azimuthal degree of freedom in the double-ring model of 
Ref. lfl9ll is irrelevant and it suffices to consider the problem 
of a BEC in a one-dimensional (ID) double-well potential, as 
in Refs. J3,[3l. Generalisation to multiple wells and different 
geometries (coupled cigars or pancakes) are straightforward. 

Different derivations of effective two-mode models have 
been presented in the literature J3, EH EUl . The goal of a 
two-mode model is generally to correctly describe the ground 
and low-lying excited states of the system. The quantity that 



3 



is obtainable from the ID model and carries unambiguous 
information about the sign of the tunnel coupling is the en- 
ergy difference Aas = Ea — E& between the antisymmetric 
state with one node and the node-less symmetric state. In the 
simplest case, the tunnel coupling k is determined from the 
single-particle linear Schrodinger equation. This approach is 
commonly used when deri ving the fully quantum-mechanical 
Bose-Hubbard model II 15l 12211 and was the basis of Ref. |7|]. 
In this case both the sign and the value of k are completely 
independent of particle number or interaction strength. The 
node theorem of quantum mechanics l23ll guarantees that the 
node-less symmetric state in a one-dimensional double-well 
potential must be the ground state, thus Aas > and conse- 
quently, the correct sign of k is negative. 

In a more general class of models based on mean-field the- 
ory, the parameters of the two-mode model are chosen in or- 
der to reproduce Aas as found from a one-dimensional GP 
equation. The ordering of eigenvalues of the GP equation by 
the number of nodes in the wave function is now no longer 
guaranteed by the node theorem of linear quantum mechan- 
ics and we are not aware of a non-linear generalization of this 
theorem. However, we find by numerical calculation that the 
ordering is preserved under repulsive interactions. The main 
result of this section is the dependence of Aas on the non- 
linear interaction strength g, shown in Fig.Q] As can be seen 
from Fig.[U the presence of a repulsive non-linear interaction 
does not change the sign of Aas and therefore k remains neg- 
ative. We now present details of our calculation. 

Starting from the three-dimensional GP equation for a BEC 
in a double-well or double-ring trap and employing a separa- 
tion ansatz, an effective ID equation describing the dynamics 
in the direction perpendicular to the potential barrier can be 
derived: 




50 100 150 200 250 300 

8 

FIG. 1: Energy difference Aas = Ea — Es between the lowest 
antisymmetric and symmetric eigenstates of a quadratic-plus-quartic 
double-well potential, plotted as a function of the dimensionless ef- 
fective interaction strength g. The fact that Aas > indicates that 
the symmetric (node-less) state remains the ground state even in the 
limit where the atoms interact strongly. Double-well parameters [see 
Eq. 0] are £o = 5 and h — 0.002 (solid curve), 0.02 (dashed curve), 
0.05 (dot-dashed curve). 



Here the energy scale Eq and length scale ao defined by the 
trap are used as units for all energies and the spatial coor- 
dinate, respectively, and the condensate wave function <f> is 
normalized to unity. We introduced the dimensionless inter- 
action strength g = gioN/ (eoao), where N denotes the num- 
ber of atoms in the trap and giu is the effective ID interaction 
strength B24ll . To be specific, we use the double- well potential 

^dw = Me 2 -£o) 2 , (7) 

where h parameterizes the barrier height between the two 
wells centered at ±£o- It is straightforward to solve Eq. (O 
with the potential (0 and find the lowest symmetric and anti- 
symmetric eigenstates as well as their respective energies Es 
and Ea- Figure [2] shows typical results obtained for low and 
high interactions strengths, respectively. As is apparent from 
the figure, the higher repulsive interaction strength is associ- 
ated with more strongly delocalized double-well wave func- 
tions, indicating an effectively stronger tunnel coupling. This 
can be explained simply by noting that the nonlinear inter- 
action energy for the two condensate fractions in each well 
shifts up their respective energies, thus effectively lowers the 
barrier and brings the condensates closer together. As a result, 
the effective tunnel coupling increases. Most importantly, the 
energy difference between the lowest symmetric and antisym- 
metric eigenstates remains positive for any strength of repul- 
sive interactions, which implies that the sign of the tunnel cou- 
pling k entering Eq. ([1]) is negative. 



IV. CONCLUSIONS 

Reference predicts a dynamical instability of a repul- 
sively interacting BEC in a double-ring trap against angular 
momentum fluctuations. We have carefully revisited the anal- 
ysis of Ref. lfl9ll and have recalculated the elementary exci- 
tation spectrum. This leads us to a different conclusion that 
makes physical sense. The ground state of a non-rotating 
condensate in the double-ring configuration is stable against 
spontaneous angular momentum oscillations. However, the 
antisymmetric state with its circular node between the two an- 
nular quantum wells can be viewed as the analog of a sta- 
tionary 2D dark soliton, which is known to have a dynamical 
instability towards the formation of local vorticity ("snake" 
instability) [20]. 
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FIG. 2: Lowest symmetric (solid curve) and antisymmetric (dashed curve) condensate wavefunctions obtained for a double-well potential [see 
Eq. 10] with £o = 5, h = 0.05, and g — 30 (left panel) or 300 (right panel). Notice the greater derealization of atoms between the two wells 
when the interaction strength is high. This arises because repulsive interactions result in an effective lowering of the tunnel barrier. 
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